On the nature of the plasma equilibrium 
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We calculate the energy of a homogeneous one component plasma and find that the energy is lower 
for correlated motions of the particles as compared to uncorrelated motion. Our starting point is the 
conserved approximately relativistic (Darwin) energy for a system of electromagnetically interacting 
particles that arises from the neglect of radiation. For the idealized model of a homogeneous one 
component plasma the energy only depends on the particle canonical momenta and the vector 
potential. The vector potential is then calculated in terms of the canonical momenta using recent 
theoretical advances and the plasma Hamiltonian is obtained. The result can be understood either 
as due to the energy lowering caused by the attraction of parallel currents or, alternatively, as due 
to the inductive inertia associated with the flow of net current. 
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Our theoretical understanding of matter is largely 
based on the Coulomb interaction between charged par- 
ticles. For small systems one can usually assume that the 
effects of the magnetic corrections are secondary but for 
larger systems this it not the case. Including magnetic 
interaction into the theories, however, has met with con- 
siderable difficulties, from the 1939 "magnetische Katas- 
trophc" of Welker [1] to the 1999 rigorous proof of the 
instability of matter with magnetic interaction by Griese- 
mer and Tix [2]. Their result is valid whether the mag- 
netic interaction is mediated by the (Darwin -) Breit po- 
tential or via a quantized radiation field. In both these 
cases it is the attraction of parallel currents that causes 
the problem. This attraction, which is so fundamental 
that it is used in the definition of the ampere, the unit of 
electric current in the Si-system, indicates a long range 
energy lowering due to correlation of currents that seems 
to diverge in many systems. 

Laboratory and astrophysical plasmas are observed to 
harbor intense currents and magnetic fields [3], but in 
plasma physics it is usually assumed that this is due to 
non-equilibrium, and that the equilibrium plasma is de- 
scribed by the traditional Maxwell-Boltzmann distribu- 
tion (see e.g. Burm [4]). This is clearly at odds with the 
above findings of an instability of the energy minimum 
ground state to parallel current generation. Alastuey and 
Appel [5] claim that inclusion of the quantized radiation 
field removes the instability problem, in direct contra- 
diction with the findings of Griesemer and Tix [2] . Most 
plasma physicists do not seem to be aware of the problem 
even if there certainly has been a fair amount of interest 
in energy minimizing states and self organization, see e.g. 
Woltjer [6], Taylor [7]. 

Here we will show that within a simple standard model, 
based on classical electrodynamics and relativistic Hamil- 
tonian mechanics, the energy of a plasma is consider- 
ably reduced when the canonical momenta are correlated 
and thus that conclusions drawn from the traditional 
non-relativistic Maxwell-Boltzmann distribution of non- 
interacting particles, or particles interacting only via a 



Debye screened Coulomb potential, cannot be trusted. 
The model, which neglects radiation, gives a quantita- 
tive estimate of the energy reduction, but does not lead 
to any un-physical divergence. 

Let us start from the following expression for the en- 
ergy of a system of classical charged particles and elec- 
tromagnetic fields, 
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(E 2 + B 2 )dV. (1) 



If we think of the particles as confined to a finite volume 
one finds that energy leaves the volume due to electro- 
magnetic radiation (Poynting's theorem). If one neglects 
the radiation as a higher order process (~ c -3 ) one finds 
that, 
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(4>Q+-A'j)dV, (2) 



is a conserved energy [8]. Inserting here that, 
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Q = £e„£(r-r (i)), j = ^e a v a (t) 8(r-r a {t)), (3) 
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one obtains, 
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for the conserved (Darwin) energy, Ed, of a closed non- 
radiating system of charged particles. Here <fi and A are 
the scalar and vector potentials, in the Coulomb gauge, 
as determined from the positions and velocities of the 
particles of the system. These can be found from the 
approximately relativistic Darwin approach [8-12]. From 
the Darwin Lagrangian, 
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where, 

L m (v a ) = -m a c 2 ^/i - vl/c 2 , (6) 

and, 



N 



K(r)= £ 



eb 



b&a) 



r - r b \ 



N 



Mr) = E 

b&a) 



and, 
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where, 

r b = (r - r b )/\r - r b \, (9) 



one also finds the expression, 



P a - —A a (r a ), 



(10) 
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where, p a — dLo/dva, is the canonical momentum. 

Since the Darwin approach normally is considered ac- 
curate to order (w/c) 2 it is sometimes stated that for 
consistency one must also expand the relativistic particle 
energy expression, E m (v a ) = m a c 2 / \/l — v 2 /c 2 , to this 
order. There are, however, several reasons not to do this. 
One is that whereas such an expansion makes the parti- 
cle energies qualitatively wrong at high speeds there is no 
reason to assume that the other terms in (4) become qual- 
itatively wrong at high speeds [12]. In fact one can show 
[13] that small modifications produce potentials that are 
valid to arbitrary speeds. A second reason is that we are 
here going for the Hamiltonian formalism and in this for- 
malism speed is not a relevant variable. The problems 
of finding the Hamiltonian corresponding to the Darwin 
Lagrangian for large systems is instead governed by a dif- 
ferent dimensionless parameter, Nr e /R, where N is the 
number of electrons, r is the classical electron radius, 
and R a typical length of the system [14-16]. 

We thus solve (10) for the velocity, 
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and insert the result in (4). We note that here it is as- 
sumed that the vector potential, A a , at particle a ex- 
cludes the contribution from that particle, as indicated 
in (8), so that self interaction is avoided. Some algebra 
then gives the result, 
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for the energy of our system of particles. This expression 
has also been derived directly by Legendre transforma- 
tion from the Darwin Lagrangian Ld (Essen [10]). 

The energy expression Ed should be compared to the 
well known energy for particles in external fields, 
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The statistical mechanics of (the non-relativistic version 
of) this Hamiltonian predicts zero magnetic response of a 
system of classical charged particles according to the, so 
called, Bohr- van Leeuwen theorem [17-19]. The expres- 
sion (12) which is valid when the field is from the particles 
themselves is different and the Bohr-van Leeuwen theo- 
rem can not be invoked to make analogous predictions 
about its properties. 

In order to find a Hamiltonian from (12) we must ex- 
press the vector potentials A a in terms of the canonical 
(generalized) momenta p a . In the Darwin approach the 
vector potential can be obtained as solution of, 
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where j t is the transverse (divergence free) current den- 
sity [20]. One can also use the ordinary current density 
(3) and impose zero divergence by a gauge transforma- 
tion afterwards [21]. In equation (14) the second time 
derivative in the d'Alembert operator of the wave equa- 
tion has been skipped since that term is ~ c~ 2 and use of 
the Coulomb gauge means that retardation is correctly 
handled to this order already [22]. In order to get the 
vector potential in terms of the p a one might then in- 
sert, 



V a = ^-^A a {Va), 
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from (11) neglecting terms of order c -2 and higher, into 
the expression (3) for the current density. For a homo- 
geneous one component plasma (with mobile particles of 
mass m and charge e) this approach transforms (14) into, 
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and [we are assuming j instead of j t in (14)], 
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essentially without approximation. Details are given in 
Essen and Nordmark [21]. 

We now assume that we, at least for long range average 
purposes, can replace the discontinuous number density 
(17) with the smoothed volume average, 

N 

™( r )=v / J2S(r-r a )dV. (19) 

This, together with the assumption of a homogeneous, 
constant density plasma, n — constant, gives, 
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r c = and j p is given by (18). We assume, as usual, 
that the homogeneous one component plasma consists of 
particles moving against a background of smeared out 
charge of the opposite sign of e, thus ensuring charge 
neutrality. The relevant solution of Eq. (20) is [21], 
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a Yukawa, or exponentially screened, vector potential. If 
we can make this momentum space vector potential di- 
vergence free we can use it in (12) to get the Hamiltonian. 
We note that it has Coulomb singularities at the particle 
positions, r = r a , but that these vanish in the Hamilto- 
nian since there self interactions are assumed excluded. 
Only the smooth part is thus of interest. 

That exponential screening should occur in a Hamilto- 
nian describing magnetic interaction was noted by Bethe 
and Frohlich [23] already in 1933. They were not aware 
of the work of Darwin, however. Neither were Bohm and 
Pines [24] who, studying collective motion of conduction 
electrons in 1951 also arrived at such a result, directly 
from a particles plus field Hamiltonian. The modern his- 
tory of the exponential screening based on the Darwin 
approach started in 1980 with Jones and Pytte [15] who 
derived it in a Fourier transformed formalism, also for a 
homogeneous one component plasma. Their derivation 
was part of a debate on the usefulness of the original 
Darwin Hamiltonian in plasma physics. 

Noting that (22) is not divergence free Essen and Nord- 
mark [21] found the gauge transformation to the correct 
Coulomb gauge divergence free expression, 
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for the vector potential from a particle with momentum 
p at the origin. Here, 

exp(x) — (1 + x) 
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and, h(x) = 2 — 3g(x). 

(24) 



We now assume that all particles of the homogeneous one 
component plasma of constant number density n have the 
same momenta p and calculate the vector potential at the 
origin by superposing contributions of the type (23) all 
over space. It is convenient to use spherical coordinates 
assuming that the momentum p is in the z-dircction. 
The contributions from the two terms of (23) separately 
diverge so it is necessary to do the angular integration 
of the second term first. Because, g{x) + h(x)/3 = 2/3, 
their sum then gives us the finite result, 
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Note firstly that without the exponential damping such 
a vector potential would diverge violently and secondly 
that the interplay between the screening length A and the 
number density h is such that both vanish from the final 
result of the integration. 

An alternative way of finding this result is by using Eq. 
(20) directly. We must then first change the left hand side 
to the divergence free electric momentum current density 
j pt . In the present case the current density is constant, 
j p (r) — (e/m)pn 1 so it is not clear what the divergence 
free version should be. Realizing that this current density 
is constant only because of volume averaging, one can use 
a result by Crisp [25] (in his Appendix A) which states 
that the volume integral of a vector field A is related to 
the volume integral of the corresponding transverse field 
M by, 

J A t (r)dV =^J A(r)dV. (26) 

When dealing with volume averages one should there- 
fore take the transverse part to be, j pt = (2/3)j p — 
(2/3)(e/m)pn. Realizing that the corresponding, A c p = 
A p t, must then also be constant (20) gives, 
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Use of the expression (21) for A 2 then again gives (25), 
A; = (2/3)(c/e)p. 

Let us return to the energy (12) which, for our case 
of a homogeneous one component plasma of particles, all 
with the same momenta p, gives, 
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per particle, if we ignore the contribution e<j) a /2 from 
the scalar potential which averages to zero (and is veloc- 
ity and momentum independent in the Coulomb gauge). 
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FIG. 1: The upper curve is the relativistic energy per particle, 
Eo(p) for free, uncorrelated, particles. The lower curve is 
the Darwin energy Ed(p), per particle, of a homogeneous 
one component plasma with a constant momentum current 
density, i. e. with correlated momenta. The energy is in units 
of mc 2 and the canonical momentum p in units of mc. 

Here we now insert the result (25) and find after simpli- 
fication, that the Darwin energy per particle is, 

1 4- 2 p, 2 ., 

E D (p) = mc 2 jggj . (29) 

V 1 + 9m 2 c 2 

One should compare this result to the usual relativistic 
energy of a free particle, 



E (p) = mc^l + ^, (30) 

which is obtained from (28) when A(p) = 0, i.e. for a 
non-magnetic plasma with uncorrelated particle motions. 
These two functions arc plotted in Fig. 1. The connection 
between velocity and momentum can be found from Eq. 
(10). For our case, p (1 — 2/3) = mv/y/l — v 2 /c 2 , so that 
p = Zmv/y/l — v 2 /c 2 . Inductive inertia thus results in 
an effective mass three times the normal one. 

The implications of the above result are far reaching for 
many aspects of astrophysical, laboratory, and metallic 
conduction electron plasmas. The latter, which are un- 
able to correlate their momenta at higher temperatures, 
due to lattice oscillations, become able to do so at lower 
temperatures. Their energy as function of canonical mo- 
mentum then shifts to the lower of the two curves in Fig. 
1. Inductive inertia thus lowers the energy of the cor- 
related electrons as originally suggested by Frcnkcl [26], 
see also Essen [27]. For astrophysical plasmas it is not 
obvious whether the nuclei or the electrons should con- 
stitute the fixed background when making the one com- 
ponent plasma approximation. The energy reduction is 



much greater if the heavy particles are considered to be 
mobile. Kulsrud [3] points out that one can explain the 
longevity of astrophysical currents as due to their large 
inductive inertia. Our result says that this inertia in fact 
reduces the energy so that it also provides a mechanism 
for the generation of these currents. 



Recall that the approximations used to get the central 
result (29) is (i) the neglect of radiation, (ii) the neglect 
of terms of order c -2 , and smaller in going from (11) 
to (15), and (iii) in the use of the smoothed density (19) 
instead of the discontinuous exact particle density. These 
approximations seem considerably less severe than those 
normally used in plasma physics. 
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